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2h ■ ABSTRACT 

We report a reanalysis of data on the measured decay rate ratio 22 Na/ 44 Ti which were origi- 
nally published by Norman et al., and interpreted as supporting the conventional hypothesis that 
nuclear decay rates are constant and not affected by outside influences. We find upon a more 
detailed analysis of both the amplitude and the phase of the Norman data that they actually 
favor the presence of an annual variation in 22 Na/ Ti, albeit weakly. Moreover, this conclusion 
holds for a broad range of parameters describing the amplitude and phase of an annual sinu- 
soidal variation in these data. The results from this and related analyses underscore the growing 
importance of phase considerations in understanding the possible influence of the Sun on nu- 
clear decays. Our conclusions with respect to the phase of the Norman data are consistent with 
independent analyses of solar neutrino data obtained at Super-Kamiokande-I and the Sudbury 
Neutrino Observatory (SNO). 

Subject headings: Astroparticle physics; Neutrinos; Nuclear reactions; Sun: particle emission 

1. Introduction by a number of groups in experiments with a 

variety of detector types and isotopes. These 
Unexplained periodic variations in measured reports, along with the observation of a change 

nuclear decay rates have been reported recently in the decay rate of 54 Mn d uring a solar flare 



Ijenkins and Fischbach ( 20091) suggest the pos- 
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behavio rs have been re p orted in the decays 
of 3 H dFalkenberal l200ll: lLobashev et all Il999t 
Veprey and Muromtsevl2012l). 32 Si. 36 C1 (lAlburger et all 



1986; I Jenkins et al.l 120091: 1 Javorsek II et allboiO; 



Sturrock et al. 2010a, 



54i\ 

60 
90 
137 



Mn (Jenkins ct al 



Cs 



201 lallbl ; I Je nkins et al.lboia) . 
201lh. 56 Mn (lEllid ll99oT . 



Co dBaurov et all l2007t IParkhomovl feOlObL ,. 
Sr (|Parkhomovll2010bllal: [S turrock et al.ll2012bh. 



2011 



( Baurov et al.ll2007t ). 152 Eu ( 



Sieeert et al 



1998), 222 Rn (and/or its dau ghters) (|Steinitz et al 



Sturrock et aH l2012al) and 2 26 Ra (and/or 



its daughters) (ISiegert et al.lll998t iJenkins et al 



20091:1 Javorsek II et al.l2010HSturrock et al.l2010 



al. 

0b, 



2011al )bl: iFischbach et al.ll2009l ). Since these flue 
tuations have been seen by groups located at 
various sites employing different detector tech- 
nologies (e.g., gas, scintillation, solid state), it 
is unlikely that they can all be attributed to 
temperature, pressure, humidity or other "en- 
vironmental" influences on the detector sys- 
tems. This conclusion is bolstered by two 
additional observat ions: A detailed analysis 
(|jenkins et al.ll2010h of possible environmental in- 
fluences on the detectors us ed in the Brookhaven 
National Laboratory (BNL) (jAlburger et al.lll986l ) 
and the Physikalisch - Techn ische Bundesanstalt 
(PTB) (jSiegert et al.l Il998l) experiments con- 
cluded that such influences were at least an 
order of magnitude too small to account for 
the observed periodic effects. Moreover, the 
changes seen over the short duratio n of a so- 
lar flare (Jenkins and Fischbach l2009h could not 
have been generated by "seasonal" influences. 
Significantly, additional periodicities of approxi- 
mately 32 d and 173 d associated with internal 
solar rotation and oscillations have been reported, 
along with the annual periodicities, in a num- 



ber of e xperiments and analyses ([S turrock et al 
2011allbl: |Parkhomo^l2^10al ISturrock et al.ll2oiobT 



cay rates by a particle or field emanating from 
the Sun. These groups have reported results from 
analyses of decay data from experiments based 
both on Earth and in spac e, and do no t find 
such an influence. Cooper ( Cooper 20091 ) ana- 
lyzed data from the radioisotope thermoelectric 
generators (RTGs) onboard the Cassini spacecraft 
enroute to Saturn and found no evidence for a de- 
pendence of the decay rate of 238 Pu (which powers 
the RTGs) on the Cassini-Sun separation. Our 
group has recently carried out a more detailed 
analysis of the RTG data, and we have shown 
that in fact there are no conflicts between the 
Cassini res ults and the previou sly reported decay 



anomalies ( Krause et al. 20121 ). Additionally, it 



was noted that even if the RTG data analyzed 
were sufficiently go od to support the claims in 
Ref. (|Cooper 2009), there would be no conflict 
in principle between the absence of an apparent 
solar influence in the a-decay of 238 Pu, and the 
observations of periodic effects in the previously 
mentioned nuclides, nearly all o f which involv e 
/3-decays (|Christmas et all Il983t iBuncherl lioioh . 
It should be noted, however, that variations in 
some a-decay r ates were found in e xperiments by 
Shnoll, et al. (|Shnoll et al.l Il998allbh . Moreover, 



even though periodic effects indicative of a solar 
influence have been seen in /3-decays, there is no 
theoretical reason why such effects should have 
been seen at the same lev el in all /3-decays. As 
we have noted elsewhere (IFischbach et al. 2009; 



Jenkins et al.ll2010l : lBuncherll2010l) . the same prop- 
erties of /3-decaying nuclei which are responsible 
for the broad range of observed half- lives (e.g., Q- 
values, nuclear and atomic wave functions, selec- 
tion rules), would likely ensure that decaying nu- 
clides would be sensitive in different degrees to an 
additional small perturbation arising from a com- 
mon solar influence. 



Shnoll et al.ll998allbi:lBaurov et al.l2007tlSturrock et al.l The preceding discussion leads naturally to the 



2010al ). Since there are no known environmental 
effects that can produce such periodic variations in 
the efficiencies of existing detectors, the implica- 
tion is that these periodicities reflect an influence 
external to the Earth environment, such as a solar 
influe nce on the decay process itself as suggested in 
Refs. (Jenkins and Fischbach 20091: 1 Jenkins et aT 



analysis of Norman et al. ( Norman et al. 20091 ) 



20091 : ISturrock et al.l 120 1 1 al ibi, bp 1 Oallbl ) 



Other groups have questioned the conclusion 
that there could be an influence on terrestrial de- 



which is the focus of the present paper. These 
authors examined data they had taken for six dif- 
ferent isotopes (five /3-decays and one a-decay), 
and found no evidence for the periodic effects re- 
ported by other authors (see above) which would 
signal a correlation between nuclear decay rates 
and Earth-Sun distance. Even though there is no 
obvious conflict between t he Norman results and 
the conclusions of Refs. ( Jenkins and Fischbach! 
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20091: I Jenkins et all |2009t iFischbach et all l2009t 
Stur rock et al Il20 1 Oallbl . I2OI 1 allbh . since none of the 



isotopes studied by Norman et al. were the same 
as those in which periodic effects have previously 
been reported, their high-quality data are useful in 
the present context since they serve to underscore 
the importance of understanding the phase of nu- 
clear decay data sets exhibiting a periodic time 
variation. This understanding will lead in turn to 
more precise limits on the sensitivities of various 
nuclides to a possible solar influence. (In what 
follows we use the word "phase" to refer to the 
calendar day on which a maximum in the count- 
ing rate is observed in a given data set. 1 January 
at 00:00 corresponds to phase 0.00, and midnight 
on 31 December corresponds to phase 1.00). 

To understand the question of phase in the 
context of the Norman analysis, we focus on the 
decay measurement ratio 22 Na/ 44 Ti which Nor- 



man et al. exhibit in their Fig. 1 ([Norman et al 



2009). These authors have graciously provided us 
with their 22 Na/ 44 Ti measurements, and we have 
reproduced their normalized residuals in Fig. 1 
after dctrending their data to take account of ex- 
ponential decay. Norman et al. compare these 
data with two hypotheses: the null hypothesis 
in which there is no annual variation, and what 
they call the "Jenkins hypothesis," which is a pro- 
posed correlation between their data and 1/R 2 , 
where R is the Earth-Sun distance. Since 1/R 2 
has a maximum at perihelion, which occurs in 
the first few days of January (varying from year 
to year between 2 January and 5 January), the 
analysis of Norman et al. builds in this phase 
at the outset even though this is not the ob- 
served phase for any of the data sets which exhibit 
annual variations. For example, the 32 Si/ 36 Cl 
phase in t he BNL experiment was 10 February 



(±3 days) (Alburger e t al.lll986tlJavorsek II et al 



20101 : ISturrock et all l2011bl ) Hence the actual 



"Jenkins hypothesis" should reflect the observed 
phase of the existing data, which can be un- 
derstood as arising from a combination of the 
annual variation of 1/R 2 and a possible North- 
South asymmetry in th e solar emission mechanism 
(|Sturrock et al.ll2011ri . The latter effect, which 
has been ascribed to the tilt of the solar axis of ro- 



tation with respect to the ecliptic ([Sturrock et al 



2011bl ). produces an additional annual variation 
with a phase of either 8 March or 8 September. 



A phase later than perihelion is consistent with 
the findings of the analyses of data taken at two 
major neutrino observatories. Analysis of data 
from Super-Kamiokande-I finds a phase shift from 
perihelion of approximately 13±17 days, which 
puts the peak solar neutrino flux at ~16 January 



(Smv and et al. fSuper-Kamiokandc Collaboration 



2004; Hosaka and et al. (Super-Kamiokande Collaboration) 



2006). Data from the Sudbury Neutrino Observa- 
tory (SNO) show a more significan t shift from per- 



ihelio n of approximately 40 days (lAharmim et al 
20051 : iRanucci and Roverei 120071 ). which puts 
the peak solar neutrino flux at approximately 
12 February. We recall that the phase of the 
BNL 32 Si / 36 Cl data i s 10 February (±3 days 



( Alburger etall Il986t IJavorsek II et all [201 



i 



Sturrock et al.ll2011bf ) and the 36 CI data taken at 
the Ohio State Research Reactor exhibits a p hase 
of 17 February (±3 davs) (| Jenkins et al.llioia ). 
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Fig. 1. — Normalized 22 Na/ 44 Ti ratio of count 
rates versus time. The solid black curve is our 
best sinusoidal fit to the data obtained from the 
Nelder-Mead (NM) optimization. With the fre- 
quency fixed at 1.00/yr, the resulting amplitude 
is 0.00034, and the phase is 0.08 (29 January). 
See text for further details. 

The impact of phase considerations on the 
analysis of Norman et al. can be understood 
by a simple example. Consider two sinusoidally 
varying functions fi(t) = Aicos (cot — <f>i) and 
fa (t) — A 2 cos (ojt — 4>2). The correlation coeffi- 
cient C(/i,/2), defined as the normalized scalar 
product of /i (t) and f% (t) , is given by 
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c{h,h) 



(MM 

l/ill/al 



cos ( 



2)- (1) 



It follows from Eq. [T] that if we identify fi (t) 
with the data of Norman et al., and /2 (t) with 
the "Jenkins hypothesis" , then the correlation co- 
efficient between /i (t) and fi it) depends on the 
presumed phase <p2 of the "Jenkins hypothesis." In 
the ensuing analysis we allow the phase and ampli- 
tude of the Norman data to be free parameters de- 
termined by the data themselves. We then demon- 
strate that in this approach the data of Norman 
et al. do in fact favor the presence of an annual 
periodicity in the 22 Na/ 44 Ti ratio, which could be 
attributable to an external influence, possibly the 
Sun. 

2. Power Spectrum Analysis 

We have analyzed the 22 Na/ 44 Ti data pro- 
vided to us by Norman, et al., using a li keli- 



hood procedure ( Caldwell and Sturrock 20051 ) de- 
veloped originally for application to the analysis 
of solar neutrino data. This procedure is similar 
to the well-known Lomb-Scargle method ( Lombl 
1976c IScargldll982l ). but it provides estimates of 
the amplitudes and phases of oscillations as well as 
the powers. In contrast to the analysis of Norman 
et al., where the null hypothesis was compared to 
an assumed sinusoidal variation whose amplitude 
and phase were fixed at the outset to be 0.0015 
and 3 January, respectively, our analysis allows 
both of these parameters to be determined by the 
data themselves. 

We formed the power spectrum of the data 
shown in Figure [TJ using this likelihood proce- 
dure. In a power spectrum analysis, the proba- 
bility of the null hypothesis (i.e., that the peak 
is due to normally distributed random noise) is 
given by e~ s , where S is the value of the power 
at a specified frequency. The resulting power spec- 
trum, over the frequency range 0-20 yr _1 , is shown 
in Figure and the strongest peaks in the fre- 
quency range 0-20 yr _1 are listed in Table 1. The 
strongest peak (with power 5=8. 35) is found at 
frequency 16.64 yr _1 . There is no obvious inter- 
pretation of this frequency, but it does hint that 
there may be unknown influences on the decay 
process. 



The peak at 1.09 yr _1 may be identified as an 
annual oscillation. We find that the amplitude of 
the normalized data at this frequency is approx- 
imately 0.00034, which is smaller by a factor of 
approximately 4 than the amplitude assumed by 
Norman et al. (0.0015). This explains why Nor- 
man et al. found a closer fit between their data 
and a flat curve than between the data and a sinu- 
soidal oscillation of amplitude 0.0015. This result 
also underscores the previous remarks to the ef- 
fect that there is no theoretical reason to expect 
that all isotopes should exhibit sinusoidal varia- 
tions of the same amplitude. We show in Figure 
the power as a function of phase, for a frequency 
of exactly 1.00 yr~ 4 . We find that the peak occurs 
at an annual phase of 0.08±0.07, corresponding to 
the date range 29 January (±24 d), i.e. to the la- 
range 5 January to 22 February. Although this is 
compatible with the phase assumed by Norman et 
al. (3 January), it allows a much broader range 
of phases, which from Eq. Q] clearly affects the 
correlation between their data and the "Jenkins 
hypothesis." We will return to this point below. 

3. Significance Estimate 

We now estimate the significance of the annual 
oscillation found by our power-spectrum analysis. 
The power at exactly 1.00 yr~ 4 is found to be 4.08, 
and we can estimate the probability that this peak 
could have a risen by chance by using the standard 
shuffle test (jBahcall and Press! 1 199 II ). We gener- 
ated 10,000 simulations of the data by randomly 
relating the actual measurements and the actual 
times. The result is shown in histogram form in 
Figure |4] We find that, of these 10,000 simula- 
tions, only 127 have powers as large or larger than 
4.08, the power of the actual data. Hence the an- 



Table 1: Most prominent frequencies in the range of 0- 
20 yr _1 obtained from the power spectr um analysis of 
the 2 2 Na/ 44 Ti data of Norman, et al. (| Norman et alj 
120091 ) 



Frequency (yr 1 ) 


Power (S) 


16.64 


8.35 


19.79 


4.59 


1.09 


4.28 


3.14 


3.24 


8.85 


2.71 
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Fig. 2. — Power spectrum of the 22 Na/ 44 Ti data of 
Norman, et al.[31]. 



nual oscillation revealed by power spectrum anal- 
ysis appears to be significant at the ~1% level. 

4. Discussion 

It follows from the preceding discussion that 
there is in the end no conflict between the results 
of our analysis and those obtained previously by 
Norman et al. The difference between our results 
and theirs is that we have addressed different ques- 
tions: Norman et al. asked whether their analysis 
of a particular dataset (the ratio of 22 Na/ 44 Ti) 
yields the same annual oscillation as might be ex- 
pected from a phase due to the 1/i? 2 variation of 
the Earth-Sun distance with an amplitude fixed 
at the outset by the BNL and PTB data. By 
contrast, we have asked whether there is evidence 
for an annual oscillation, for any combinations of 
amplitude and phase. As shown in our Figure 
there is a wide range of values for the amplitude 
and phase for which an annually varying signal 
provides a better description of the Norman data 
than does the null hypothesis. 

We therefore adopt the position that since the 
mechanism of these oscillations is currently un- 
known, and since our analyses of several datasets 
yield evidence of annual oscillations with a range 
of amplitudes and phases, there is no reason to 
expect that the amplitude and phase of an annual 
oscillation in the Norman dataset should be iden- 
tical to that of the amplitude and phase of either 



Fig. 3. — Power versus phase for annual oscilla- 
tion. 

the PTB 226 Ra measurements, or to the ratio of 
the BNL 32 Si and 36 C1 measurements. This ob- 
servation can be strengthened by noting that al- 
though one can identify an annual oscillation in 
each of the BNL 32 Si and 36 C1 measurements, the 
individual phases of each nuclide differ from the 
phase of their ratio and each other: The phase 
of the 32 Si measurements is 0.899 (24 Nov) and 
the p hase of the 36 CI meas urements is 0.695 (II 
Sep) (jSturrock et al. 2011bl ). This fact indicates 
that the phase is probably not determined solely 
by the Earths orbital motion and emphasizes the 
importance of leaving the phase as a variable to be 
determined from the data, for each set of measure- 
ments. We also note in passing that since the BNL 
32 Si and 36 C1 data were acquired at the same time 
on the same apparatus, the very fact that these 
data exhibit different phases argues against at- 
tributing the observed annual variations in count 
rates to an environmenta l influence on the appa- 

l20inh . 



ratus (Jenkinsetal 



Our focus here as been on annual oscillations 
in some measured nuclear decay data, since these 
were the only variations discussed by Norman et 
al. However, these are not the only periodic- 
ities which have been observed, and they may 
not be the most signif i cant. In recent articl es 
dJavorsek II et all l2010t JSturrock et al.l l2011al lbl. 
2010bUal : Fischbach et al.ll2011 ) , we have presented 
evidence for periodicities in nuclear decays with 
frequencies in the range 11-13 yr _1 and ~2.11 




x 10 



Fig. 4. — Histogram of the power at 1.00 yr _1 , 
formed from 10,000 simulations, generated by the 
shuffle procedure. Only 127 of these simulations 
have power as large or larger than the power 4.08 
of the actual data shown by the vertical line. 



yr . The 11-13 yr _1 periodicities may be at- 
tributable to internal rotation of the Sun. The 
~2.11 yr -1 periodicities may be an analog of the 



well-known Rieger oscillation (IRieger et al 



19841) 



attributable to an r-mode oscillation (|Saiol 11982) 
in an inner solar tachocline. It would clearly be of 
interest to extend this and other analyses of decay 
data to look for such periodicities which, as noted 
above, cannot be attributed to known environmen- 
tal effects. Furthermore, an examination of the 
literature yields additional evidence for these ob- 
served periodicities, annual and otherwise, which 
have also been found in decay data. A list of the 
experiments we have found to date is presented in 
Tabled 

As we show in Table [2] the majority of the 
isotopes listed have at least one beta-decaying 
branch, or several beta-decaying daughters. There 
is one exception, however, an d that is the re- 
sult presented by IShnoll et al. I (|l998alfbh . where 



an oscillation was observed in 239 Pu, which is 
strictly alpha-decay through any realizable daugh- 
ter that could be measured. However, others have 
made measurements of isotopes which decay by 
alpha with no realizable beta-emitting daughters 
a nd have seen no such effe cts, e.g. th e result s 
of Norman et al. | (120091) and lParkhoinovl (l2010bl) . 
The mechanisms of beta- and alpha-decays are ob- 
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Fig. 5. — Plot of the difference in reduced chi- 
squared values (null hypothesis "Jenkins" ) with / 
= 1.00 yr -1 , for various amplitudes and phases. 
Colored regions correspond to sinusoidal parame- 
ters where the "Jenkins hypothesis" is preferred. 
The optimal ht, located in the central, darkest re- 
gion with amplitude « 3.4xl0 -4 and phase ~29 
d, has a reduced chi-squared difference of «0.14. 
1 January at 00:00 corresponds to phase 0.00, and 
midnight on 31 December corresponds to phase 
1.00. 

viously very different. While we still have not de- 
termined a mechanism by which solar neutrinos 
would affect the weak interaction associated with 
beta-decays, the development of a model where 
the neutrinos could affect both alpha- and beta- 
decays becomes even more difficult. We note, how- 
ever, that there are thousands of isotopes of vary- 
ing decay modes which have not been examined. 
Many of these measurements will be complicated 
due to very short half-lives, or other properties, 
which may make many measurements impractical 
with present technologies. This includes the beta- 
and gamma-decays of all known elementary par- 
ticles, including the neutron, and many nuclides. 
However, there remain many avenues to explore in 
this work. 

We are indebted to E.B. Norman and his col- 
laborators for generously sharing with us the data 
analyzed in this article. The work of P.A.S. was 
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06072572, and that of E.F. was supported in part 
by U.S. DOE Contract No. DE-AC02-76ER07128. 
The views expressed in this paper are those of the 
authors and do not reflect the official policy or po- 
sition of the U.S. Air Force, U.S. Department of 
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Defense, or the U.S. Government. The material in 
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Force Flight Test Center Public Affairs. 
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Table 2 

Some experiments where time-dependent decay rates have been observed. 





Decay 


Detector 


Radiation 


Isotope 


Type 


Type 


Measured 


3 H 


P- 


Photodiodes 


P~ 


3 H 


P~ 


Liq. Scint. 


P~ 


3 H 


P~ 


Liq. Scint. 


P~ 


3 H 


P- 


Sol. St. (Si) 


P~ 


22 Na/ 44 Ti a 


/3+,K 


Solid State (Ge) 


7 


36 CI 


P- 


Proportional 


P~ 


36 CI 


P~ 


Geiger-Miiller 


P~ 


54 Mn 


K 


Scint. 


7 


54 Mn 


K 


Scint. 


7 


56 Mn 


P~ 


Scint. 


7 


60 Co 


P~ 


Geiger-Miiller 


P-,1 


60 Co 


P~ 


Scint. 


1 


85 Kr 


P~ 


Ion Chamber 


7 


90g r/ 90 Y 


P~ 


Geiger-Miiller 


P~ 


108m^g 


K 


Ion Chamber 


7 


133 Ba 


P~ 


Ion Chamber 


7 


137 Cs 


P- 


Scint. 


7 


152 Eu 


P-,K 


Sol. St. (Ge) 


7 b 


152 Eu 


P~,K 


Ion Chamber 


7 


154 Eu 


P~,K 


Ion Chamber 


7 


222 Rn c 


a,p- 


Scint. 


7 


226 Ra c 


a,P~ 


Ion Chamber 


7 


239p u 


P~ 


Sol. St. 


a 



Reference 



Falkcnbcrg (2001) 

Shnoll et al. (1998a .b") 

Vcprcv and Muromtscv (2012) 

Lobashev et al. (19991 

Norman_et_aL (2009) and this article 

Jenkins et al. (2009) ; Sturrock et al. (2010a . 2011a' 

Jenkins et al. (2012) ™ 

Jenkins and Fischbach (2009) 

(Jenkins et al. 2011 ) 

Ellis (1990) 

Parkhomov (2010b .a) 

Baurov et al. (2007) 

Schrader (2010) 

Parkhomov (2010b .a): Sturrock et al. (2012b) 

Schrader f2010) 

Jenkins, et al., arXiv 

Baurov et al. (2007) 

Siegert et al. (1998) 

Schrader (2010) 

Schrader (2010) 

Steinitz et al. (2011) ; Sturrock et al. (2012a) 
Jenkins et al. (2009) ; Sturrock et al. (2010b . 2011a 
Shnoll et al. fl998a .b) 



a Only the count rate ratio data were available. 
b Only the k photon was measured. 

c Decay chain includes several primarily /3-decaying daughters which also emit photons. 
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